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The leading order low energy parameters like the pion decay constant or the quark condensate 
are well-known from "classical" low energy theorems and experiments. It is a challenge, how- 
ever, to find these parameters based exclusively on an ab-initio QCD calculation. We discuss 
results of a quenched lattice calculation of low energy constants using the chirally improved 
Dirac operator. Several lattice sizes at different lattice spacings are studied, using pseudoscalar 
and axial vector correlators. We find consistent results for f K = 96(2) MeV, /k = 105(2) MeV, 
£=-(286(4) MeV) 3 , the average light quark mass m = 4.1 (2.4) MeV and m s = 101 (8) MeV. 
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1. Setup 

When QCD was accepted as the quantum field theory of strong intereractions it became clear 
that its (approximate) chiral flavor symmetry is spontaneously broken. The exploitation of the 
underlying principles has then led to the development of chiral perturbation theory. In that context 
a systematic expansion of many observables in terms of low energy constants has been derived. 
These constants, however, have to be determined either from experiment or from basic principles, 
i.e., non-perturbative solution of the underlying field theory QCD. The lattice formulation of QCD 
allows such a determination. 

Since chiral symmetry and its spontaneous and explicit breaking play an essential role, it is 
desirable to work with fermionic actions having that symmetry. Lattice Dirac operators obeying 
the Ginsparg-Wilson relation (GWR) are the lattice analogue of chirally symmetric continuum 
Dirac operators. Here we use the chirally improved Dirac operator Dq\ [jl|]. It is based on a 
systematic expansion of the lattice Dirac operator taking into account the whole Clifford algebra 
and terms coupling fermions within a certain range of neighbors on the lattice. The expanded Dirac 
operator is inserted in the GWR which then leads to a set of algebraic equations for the expansion 
coefficients. In recent applications [fj], as well as here, we use a set of 19 independent terms in the 
action. Finally, in the definition of Dqi also one step of HYP-smearing of the gauge configuration 
is included. For the gauge fields we use the Liischer-Weisz action. The lattice spacings have been 
determined using the Sommer parameter. We summarize the simulation parameters in Table [T[ 

In Ref. [Q] hadron masses for a quenched simulation based on Dq\ have been presented. Pre- 
liminary results for low energy constants in that context have been published in Ref. At that 
time the renormalization constants relating these to the continuum MS -scheme were not available. 
Meanwhile the necessary constants for quark bilinears have been determined for Dei in [Q]. This 
now allows us to compute some of the basic low energy parameters in the quenched case. For this 
purpose we study quenched QCD at various values of the quark masses and determine results for 
m u = ma ^ m s in the ud and us meson sector for several lattice sizes and lattice spacings, down to 
a pion mass of 330 MeV. 

All our results have been derived from correlation functions of pseudoscalar interpolating 
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Table 1: Parameters of the simulation. Where the strange quark mass is given (in parenthesis), we also 
determined propagators for sttange hadrons. Type denotes the type of the quark source/sink and #cf. the 
number of configurations entering the analysis. 
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Figure 1: my> vs. Mj, s using X = P. Here Mp$ denotes M K or M%, respectively. All masses are given in 
physical units in the MS-scheme. 

fields P = dy=,u and A4 = dy^y^u. Where we give physical, renormalized (MS-scheme) values, 
we always use the renormalization factors as obtained in the chiral limit, even in plots for non- 
vanishing bare quark masses. The renormalization factors Z relate to the continuum MS-scheme at 
a scale of \i = 2 GeV and are a-dependent as given in [|], §]. We use the interpolating fields P and 
A4 for both the pion (with degenerate light quark mass rn), and the kaon, with one light and one 
strange quark. The bare mass of the heavy quark is fixed with the help of the physical kaon mass. 
In order to improve signals we use Jacobi-smeared sources and sinks for the quarks. 

2. Results 

In Ref. [§] we discuss the normalization of the sources, methods for fitting the propagators, 
error estimates (based on a jack-knife analysis), numerical derivatives, and the results for the meson 
masses. Here we present only some results for quark masses, condensate and the pseudoscalar 
decay constants. 

2.1 Quark masses 

The renormalized quark mass may be determined, utilizing the axial Ward identity, 

Z A (d t A 4 (p = 0,t)X(0)) Mlf K 



Z P (P(p = d,t)X(0)) 



■ Z m 2m = 2m (r> , (2.1) 



where X may be P or A4. The asymptotic behavior of the correlators cancels in this ratio and 
the plateau values provide the mass (proportional to the so-called AWI-mass). For the second 
choice, X = A4, the correlators are of sinh-type and the ratio becomes numerically unstable near 
the symmetry point in t. Our analysis of the quark masses is based on the more stable choice X = P. 

The quark mass data in Fig. |l|, presented in physical units in the MS-scheme, leads to a very 
consistent picture. The abscissa gives the corresponding pseudoscalar mass, i.e., that of the pion or 
the kaon, for the corresponding values of mV'. We find that for given lattice spacing the numbers 
for kaon and pion are on top of each other, although these states have different quark content. 
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Figure 2: Scaling properties of the light quark condensate (l.h.s.) and of the pion decay constant (r.h.s). 

The lines in Fig. [j] correspond to linear fits mS r ' °< M\, enforcing the simultaneous chiral limit 
of both observables. The data do not show deviation from that behavior, although logarithmic 
corrections are expected due to quenching. The linear extrapolations are in good agreement with 
the Particle Data Group [|6|] average for the light quark and strange masses at the physical pion 
mass. We obtain 

^ (m[ r) + m { f} = fflW ~ 4. 1 (2.4) MeV , ^ (m { s r) + m w ) -52(3) MeV (2.2) 

in the MS-scheme. Combing theses values gives m| r ' = 101(8) MeV. Possible finite size effects 
and other systematic effects like chiral extrapolation and quenching have not been accounted for. 
The given error takes into account the standard error and includes the derivations due to the de- 
pendence on the lattice spacing. The numbers are in good agreement with determinations from the 
overlap action in [f7], §]. 

2.2 Chiral condensate 

We have computed the renormalized condensate from the relations 

j*At* = -2mW>, (2.3) 
Z A Z P (A A (p = 0j)P(0)) f^*' = ^ e - M «< = |EW| e -^, (2.4) 

Z A Z P yj (A 4 (p = 6,f )A 4 (0)> (P(p = 0,t)P(0)} ~ |EW| e -^«' , (2.5) 

which all contain EM. The first of these is the GMOR relation, the other two are determinations 
directly from the coefficients of propagators and implicitly related to GMOR as well. 

We find excellent agreement for all three determinations; the values are consistent within the 
error bars. The dependence on the bare quark mass is compatible with the leading (linear) chiral 
behavior. Note, that we are not at small enough quark masses to be in the so-called £-regime [^] 
but are in the p-regime. 

In Fig. ^| (l.h.s) we show the results of the linear extrapolation to the chiral limit for all three 
types of determination for all lattice sizes studied. The expected scaling behavior &(a 3 ) is demon- 
strated by the fitted curves. Assuming the leading scaling behavior throughout (as was done in 
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Figure 3: af K x vs. (aMps) 2 with error bands for chiral extrapolations as discussed in the text. The fit 
includes only data in the indicated range. 

other studies where only one lattice spacing was studied) the average of the resulting values for 
the condensate in the continuum limit is |EW| = (286(4) MeV) 3 . This is slightly larger than a 
determination from the overlap action |jl2|, [7|] (although still with in the error limits) and larger than 
the corresponding results in Ref. [[| 13]. Our numbers are in good agreement with calculations in 



the e-regime 0, [TT]]. 
2.3 Decay Constants 

The pseudoscalar decay constants have been extracted from the asymptotic behavior of the 
pseudoscalar correlation functions according to 

Z\ (A 4 (p = 0,t)A 4 (0)) r,Mpsf^- Mpsl (2-6) 

for pion and kaon, respectively. When plotting them as functions of the respective pseudoscalar 
masses in Fig. Bl the data for pion and kaon essentially overlap each other and exhibit a universal 
functional behavior. We also show the error band of a quadratic extrapolation to the chiral limit. 
In quenched QCD one expects correction terms with a logarithmic singularity in the valence quark 
mass m. As pointed out in Jl4]], the leading order logarithmic term m logwi of ChPT involves quark 
loops that are absent in the quenched case. There will be non-leading, e.g., logarithmic, terms, 
though. In addition to the term linear in the quark mass m we therefore also allow a term m 2 log m 
in the extrapolating fit (cf. the discussion in [|I5|]). Actually, in the fit it makes no significant 
difference whether we take this term or just m 2 . 

In Fig. ^|(r.h.s.) we show the scaling behavior by comparing results for different lattice con- 
stants for our largest physical size lattices. The results are from the chiral extrapolation for the 
light quarks; in the (semi-) chiral extrapolation for the kaon decay constant the strange quark mass 
parameter is held fixed. Assuming the leading @{a) behavior, the continuum limit values are 

fn = 96(2) MeV , f K = 105(2) MeV . (2.7) 



Studies for the overlap action in quenched simulations have obtained similar results for f n [JT5j, [T 
IIJ 0@ [Dp. The experimental values are /„ = 92.4(0.3) MeV and f K = 113.0(1.3) MeV [§], 
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